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Abstract 

We have demonstrated that it is possible to access a crossover scenario starting with a weak 
couphng (BCS) d-wave superconductor to a strongly coupled Bose-Einstein condensate (EEC) 
phase as the exchange interaction is tuned in a two dimensional system described by a t-J-U 
model via numerically solving the Bogoliubov-de Gennes (BdG) equations. While in the extreme 
dilute limit, the electronic pairing phenomena is independent of the Coulomb repulsion, U, the 
superconductivity depends on U, and so does the crossover. Further, the effect of variation in the 
carrier density on the BCS-BEC crossover has also been investigated. The crossover picture is 
illustrated by computing the chemical potential, which when falls below the noninteracting band 
minimum, signals the onset of a phase with tightly bound, shorter pairs. As an evidence of the 
above feature, the Cooper pair radius is calculated which shows a significant shortening at the 
emergence of a BEC-like phase. Besides, in contrast to the previous work where a crossover was 
claimed only in the dilute limit, we have demonstrated it at large densities near half filling. 
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I. INTRODUCTION 



The phenomena of BCS superconductivity to Bose-Einstein condensation crossover is an 
old one[l|-l3[ and is interestingdue to its possible connection to pseudogap observed in high- 
Tc superconductors (HTSC) [4-^ and also understanding the rich phase diagram of these 
materials, j?, 8| Experimental data suggest that the pairing symmetry in the HTSC cuprates 
is predominantly of dx2-y2 wave symmetry, jol. [lO | Unlike the conventional BCS superconduc- 
tors, these compounds are characterized by extremely small pairs having a spatial extension 
of the order of one (or a few) lattice spacing(s). Randeria and coworkers shown that in two 
dimensional (2D) systems, a two-body bound state in vacuum is necessary for an s-wave 
(/ = 0) pairing instability, which is not true for the higher angular momentum channels 
(/ > 0). ll| The threshold for a single bound state of two electrons on the empty lattice of 
a 2D t-J model are = 2t and = 2t/(4/7r — 1) ^ 7.32t for s- and (i-wave symmetry, 



respectively. 
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The BCS superconductivity emphasizes large, overlapping Cooper pairs formed by weak 
attractive interaction among the electrons. With increasing the attractive potential, the 
Cooper pairs smoothly evolve into nonoverlapping composite bosons, the so called BEC 



with larger binding energy. 1^, [l5| Earlier studies mainly focused on the evolution of the 

-[l9| The combined effects 
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coherence length of a pair of fermions at zero temperature, 
of density and interparticle potential show that crossover is robust for all densities in case 
of s-wave pairing and only for low densities in case of rf-wave pairing. 23|] Further, in 2D 
(i-wave systems, the existence of a finite range of potential and the next-nearest-neighbor 



3, 



25l | Experimentally measurable ther- 



3-22] Ahhough, the igated from 

various different aspects, however a conclusive study on the BCS-BEC crossover in a d-wave 
SC is still pending. 

Coulomb repulsion enforces the no-double-occupancy constraint and has been proposed 
to be the key ingredient for HTSC. 26|] Recently the t-J-U {U > 0) model, known as 



hopping drastically influence the crossover diagram. 

modynamical properties have also been calculated as a function of density and interaction 
strength across the crossover. 



the Gossamer Hami 
diagrams of HTSC. 



tonian has been successfully used to describe the ground state phase 



27H33| At half-filling, the Gossamer SC state undergoes a quantum 



phase transition to an antiferromagnetic Mott insulator with increasing on-site Coulomb 
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repulsion. 
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lence bond SC state 



3l| Near half-filling, the Mott insulator evolves into a resonating va- 



32| as we increase the Coulomb repulsion. Further, a phase separation 



between antiferromagnetic and rf-wave superfiuid phases appears for U > 7.3t, 33] which 
indicate the role of Coulomb repulsion on spin and charge density waves and stripes in 
cuprates. Apart from describing the competing orders, t-J-U model is useful for a better 
understanding of electronic inhomogeneity and local electron density of states observed in 



scanning tunneling microscopy measurements. 



34l ] f/ — 7- oo limit projects out all the doubly 



occupied sites and renders a more familiar variant of the model, normaly the t-J model. 

Considering the important role of on-site Coulomb repulsion in HTSC, we investigate 
its impact on the critical pairing strength for bound state formation and on the BCS-BEC 
crossover near half filling. We use the t-J-U model and recalculate the threshold for two 
particle bound state ~ 7.32t for completeness. |l2[] This result is independent of U as it 
should be, since U has no role in a two particle intersite pairing. In fact, this result will be 
contrasted for pairing in a many-body system. For the many-particle case, we employ the 
self-consistent Bogoliubov-de Gennes (BdG) method to compute the relevant physics. We 
define Jsc as the threshold below which rf-wave correlations are zero, and the significance of 
which is illustrated later. We find that Jg^ initially increases with density and then decreases 
closer to half filling (n=l), showing a characteristic optimization behavior. We also find the 
signature of BCS-BEC crossover close to half filling and show that the chemical potential 
crosses the noninteracting band minimum with increasing the pair correlation strength. Our 
claim is well supported by the results obtained for the Cooper pair radius. It is worth while 
to point out that while the BdG studies on s-wave superconductors are fairly abundant, 
similar attempts on (i-wave systems are limited. The problem is comprised by the fact that 
the attention that rf-wave systems recieved is mostly restricted by very narrow regime of the 
parameter space. We believe that attainment of the self-consistency is the most crucial issue 
in all these studies, owing to, possibly a large number of competing ground states. We have 
been able to address this difficulty partially by a very thorough numerical investigation of 
the problem. 

Organization of the paper is as follows: a brief introduction to the model Hamiltonian is 
given in Section-II. In Section-Ill, we calculate the critical exchange potential, Jgc for two 
particles in 2D using the t-J-U model. In Section-IV, we consider finite particle density and 
investigate superconducting pairing in general, and the impact of on-site Coulomb repulsion 
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on the pairing correlations, and finally on BCS-BEC crossover by calculating the chemical 
potential and mean-pair radius. Finally, we conclude in Section-V. 



II. MODEL HAMILTONIAN 



We model the 2D (i-wave SC by the Hamiltonian, 

H = HkE + 'Hint + 



(«i>:0 (ij) 

+t/ ^ ni^n^ - /i ^ ni . (1) 

i i 

The first term describes the kinetic energy, where t is the hopping amplitude of electrons 
between nearest-neighbors (ij) with spin a, J is the pairing interaction strength between 
the nearest neighbour sites and U is the on-site repulsive interaction. It is easy to show that 



the interaction terms lead to a (i-wave SC ground state in the singlet channel. [12| The spin 
operators are defined by Sf = cl^o''^i^Cii3, where the cr° are Pauli matrices, and the density 
operators Uia = cj^Cia with rii = rii^ + rijj^. The chemical potential /i is adjusted to fix the 
average carrier density n = 



III. TWO PARTICLE BOUND STATE 



For the sake of completeness and to compare and contrast with the results for finite 
electron density, we present the calculation of threshold exchange interaction in the extreme 
dilute limit, that is for two electrons in an empty lattice. This specifies the condition for a 
two-particle d-wave bound state. 

Here we start with the particle wave function that is appropriate for a singlet pairing, 

|v[/) = 5^$(^i,^2)4^4j0) , (2) 

where $(21,22) = ^{i2^H)- Thus the equation of motion, 'H|\E') = £\^) can be written as {£ 
being the energy of the electron pairs), 

8^{iui2) = ^[t,,j$(j;z2) +ti2i$(«i,j)] + [f/5n,i2 - ^n,ij$(«i,^2) (3) 

4 



Fourier transform of the equation yields, 



U 



£<l>(ki,k2) = [t(ki)+t(k2)]$(ki,k2) + -5^$(ki+k,k2-k) 



where 



$1 



-l5^J(k)<|.(ki-k,k2+k) , 

k 

ki,k2) = j^^'^ih^k) X exp[-i{ki ■ r;, + ka ■ J] 



(4) 



(5) 



t(k) = — 2t(cosA;a; + cos ky) and J(k) = 2J{coskx + cos ky), where the lattice constant is 
chosen to be unity. Let us define Q = ki + k2, q = |(ki — k2), and $(ki, k2) = $Q(q) then 
we obtain 

''E.<fQ(k)-iE.^(q-k)<fQ(k) 



*Q(q) 



N 



Q 



(6) 



E-ti^ + q)-t[f-q 
For a singlet pairing, we can take Q = and Eq. can be decoupled so that we can write 

Co = U IqCq — 2JIxCx — 2J lyCy 
Cx = UCqIx — 2JCxIxX ~ '^JCylxy 
Cy = UC^Iy — 2JCxIxy ~ '^JCylyy 

C's and J's are lattice integrals which can be defined as. 



(7) 



Co 
h 
h 



q 

N ^ 

N ^ 



Ci 



-E 



cos /c/$n(k) 



1 



^ E + 4t(cos Qx + cos Qy) ' + 4t(cos Qx + cos Qy) 



cos Qx COS Qy 



COS qi 



N ^ E + 4t(cos Qx + cos Qy) 



cos^ qi 



N ^ E + Atcosqi 



(8) 



where / G {x,y}. For an isotropic square lattice symmetry ly = Ix, lyy = Ixx and Cy = Cx- 
Unique solutions of Cq, Cx and Cy can be obtained iff the determinant of the coefficient 
matrix is zero i.e, 

UIo - 1 -2JIx -2JIx 



UIx 



-2 J Ixx - 1 



-2JL 



xy 



2 J Ixy 2 J Ixx f 







(9) 
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To project out the possibility of double occupancy we now consider U — )■ oo. This reduces 
the Hamiltonian in Eq. ([1]) to a more familiar variant, the vastly studied t-J model. It is 
worth mentioning that the lattice integrals can be evaluated analytically and one can verify 
that only two integrals are completely independent, such that we only need 

T ^ a ' 

-'0 — o 

_ l {a'-2a + 2)K{=^) 1 (a - 2)n(f , ^) , l aEj-^) 
~ 4 Tra 2 n + 4 vr ' ^ ^ 

where a is a dimensionless energy, j^. K, E and 11 denotes the complete elliptic integral of 

first, second and third kind respectively. The rest of the integrals can be represented using 

these two in the following way, 

E \ E 

hx + hy = --rr^x, and 4 = — - — -^o • (H) 

4t ot ot 

So the lattice integrals can be computed in terms of complete elliptic integral of the first 
dnd K{—2/a), the second kind E{—2/a), and the third kind Il{2/a,—2/a) respectively 



35|. The critical value of J can be obtained by substituting S = —8t — 6, i.e the enegy for 
the electron-pair lying just below the noninteracting two electron band in 2D (= — 8t) from 
Eqs. (ITOl [TTjl . Expansion of the elliptic integrals for 5—7-0 yields a logarithmic divergence. 
However the coefficient of the diverging term should be made to vanish, which leads to a 
quadratic equation for J, viz. 

(4-7r)j2 -8J + 47r = . (12) 

The solutions of Eq. ( fT2|) are as Jc = 2t and 7.32t, which actually leads to the critical value 

of bound state formation for two electrons corresponding to s-wave and (i-wave symmetries. 

Conversely, from Eq. ([9]) we obtain a pair of equations as, 

l-[//o 



J' 



2(2f//2 + (4, + 4,)(l-f//o)) 



The first equation in Eq. ( fT3l) corresponds to the s-wave case and the second one to the 
d-wave. One should note that the d-wave bound state is independent of U . The s- and 
(i-wave bound state energies can be written as, 

ttJ 



El = -8t-64texp 
Ef^ = -8t- 



J-J- 
16t2 r 1 1 



1-^ 



(14) 
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FIG. 1: (color online) Behavior of binding energy Eh/t as function of pairing strength J/t for a 
isotropic square lattice with s-wave and d-wave symmetry. The critical pairing strengths for bound 
state formation are 2t and 7.32t in case of s- and d-wave, respectively. 

where J^'*^ being the critical value of J in s- and rf-wave cases. In Fig. ([T]), we show the 
dependence of the bound state energies for the s- and (i-wave cases as a function of the 
exchange interaction, J. It can be seen that the s-wave critical value drops exponentially at 
lower J values, whereas the (i-wave curve has a near-linear fall off. However a s-wave bound 
state is energetically more favorable compared to the other. It is interesting to note that, 

is higher that J*, despite the fact the latter is affected by the strong on-site repulsion. 
Thus we conclude that for J > aX low density, the ground state of the system is a Bose 
condensate of fermion pairs having a s-wave pairing symmetry. 

Hence we analyze pairing correlation in moderately high particle density regime in the 
coming section, where a crossover to a BEC state is further investigated. Competition 
between the two interactions results in (a) a rich behavior of Jsc with increasing n and U, 
and (b) BCS-BEC crossover close to half filling, in contrast to the existing results. 
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IV. BCS-BEC CROSSOVER IN A d-WAVE SUPERCONDUCTOR 



To solve the many-particle Hamiltonian (Eq. [T]) we reduce the quartic terms to quadratic 
terms with a mean-field decoupling and then employ Bogoliubov transformation 

Cia = {^inalna " OlV-^alla) (15) 

n 

to diagonalize the mean-field Hamiltonian, which takes the form 

H = J2EnalLlna ■ (16) 

na 

7^ ^ creats an elementary fermionic Bogoliubov quasiparticle excitation with quantum num- 
ber n, spin a, and energy En^ > 0. Calculation of the commutators of the above H with 
the electron operators Qo- leads to the BdG equations 

where iun{i) = -Y.si't + Wi)un{i + 6) - fiiUn{i) and Aun{i) = Y^s^ihi + S)un{i + ^), 
and similarly for Vn{i)- The pairing amplitudes on a bond is defined by A(z; i + 6) = 
— J{ci+siCi^ + C4Q+5^)/2, where 6 = ±x,±y. The Hartee-Fock shifts are given by fli = 
11 — U (rii) + ^ '^^{nij^s) and Wi = ^(c|_,_^ _^Cj^cf)- Starting with an initial guess for all local 
variables defined on the sites and bonds of the lattice, we numerically solve for the BdG 
eigenvalues for En > and the corresponding eigenvectors The pairing amplitude 

is given by 

n 

with the density (nj) = 2^^|f„(i)p, and the Fock shift Wi = f ^„'yn(^ + S)vl{i) are 
calculated iteratively until self-consistency is achieved. The (i-wave pairing amplitude is 
defined as, 

= [A{t; +x) - A{i; +y) + A{t; -x) - A{t; -y)] /4 , 

where x and y connect to the nearest neighbors in 2D. 

It is worthwhile to mention at this point that achieving self-consistency in our calculations 
is a significantly crucial issue, as a number of parameters demand simultaneous convergence, 
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FIG. 2: (color online) Variation of critical value of pairing strength Jgc/t with density for fixed 
values of [/ (= 1 and 2). Jgc shows an optimization behavior with density due to the on-site 
Coulomb repulsion. The system size considered is 24 x 24. 



which as experienced by us, is a challenging task possibly because of a number of competing 



ground states. May be due to this difficu' 
by taking only one set of parameter values. 



ty Ghosal et. al. have done their calculations 



36| Thus an extensive search for self-consistent 



solution over a broad parameter space has been done, however with limited success. Yet in 
our work we were able to scan a reasonable amount of parameter space defined by U , J and 
n, which has not been done earlier. 

We begin by studying the ground state properties and calculate the critical value of 
pairing energy J^c which we defined as a threshold value below which pairing correlations are 
zero (normal state). We find that Jg^ depends strongly on the density and on-site Coulomb 
repulsion. Then, we investigate the impact of Coulomb repulsion on the BCS-BEC crossover 
picture. The density range we are interested in is between quarter and half filling as the 
situation for lower densities is addressed in earlier studies. We consider a 24 x 24 system 
size for all our calcuations. 

We have seen in the last section that the threshold for a single bound state of two electrons 
on the empty lattice of a 2D t-J model is ^ 7.32t for rf-wave symmetry. With increasing 
particle density the average spacing between the particles decreases and consequently the 
overlap of the single particle states increases. Thus, the minimum pairing potential required 
for bound state formation monotonically decreases with density. |23|] In a t-J-U model, as 
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FIG. 3: (color online) Behavior of d-wave correlations (a) with pairing energy J shows a linear 
grow and (b) with electron density n saturates near half filling (n=l) due to complete overlap of 
the electronic states. The parameters are taken from Fig. [2] and chosen so that the system always 
in BCS or BEG state. 

we increase density, the on-site Coulomb repulsion becomes more effective, which makes 
it energetically costlier for the electrons to come closer and form pairs. So in our case 
Jsc initially increases with density, with its value always less than 7.32t (Fig. |2]). Close 
to half filling, the electrons naturally sit next to each other resulting in lower minimum 
pairing energy. Hence the critical coupling shows a characteristic optimization behavior 
with particle density. This non-monotonic behavior stimulate our interest to investigate the 
BCS-BEC crossover near half filling, where it always shows characteristics of a BCS state 
in the absence of \J . It may be noted that even though the two-body bound is insensitive 
to the value of the Coulomb repulsion, the pairing correlation in a system with finite carrier 



density, this is no longer true. Further, as our main 



'ocus is the larger density regime, so we 



do not address densities lower than qurter filling. 23 1 

Both the attractive pairing potential and electron density drive (i-wave correlations, A^. 
As Arf is direcly proportional to the pairing potential, it grows linearly with the strength of 
the potential for a fixed particle density (Fig. |3]^a)). Further, the (i-wave correlations depend 
upon the overlap of the electronic states, so increases with particle density. Thus, for 
fixed values of J and V the correlations monotonically increase with density. Close to half 
filling, the particles start sitting next to each other and consequently the overlap becomes 
complete, resulting in a saturation behaviour of A^ with density (Fig. [^b)). 

The BCS-BEC crossover scenario is conveniently investigated by calculating the chemical 
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FIG. 4: (color online) Variation of (scaled) chemical potential, /i' as a function of (a) pairing 
energy J and (b) electron density n. BCS-BEC crossover occures near half filling {n=l) for U=2 
and moderate pairing strength. 

potential as a function of the interaction strength between the fermions^, [s^]- In general, 
bosonic degrees of freedom is expected to emerge once the chemical potential of the many- 
body ground state slips below the noninteracting band minimum in a tight-binding system. 
In a s-wave system, a smooth crossover from fermionic superconductivity to bosonic degree 
of freedom can occur for all densities as the coupling strength is increased. 23| In a (i-wave 
system, without the on-site repulsion bosonic degrees of freedom can only emerge in the 
extreme dilute limit, while for the large densities, the system behaves more like a weak- 



coupling superconductor 



231]. Although, as mentioned erlier, the on-site Coulomb repulsion 



is not indispensably necessary for generating a two-body rf-wave pairing state, we find that 
it has a significant role in (i-wave pair formation and the BCS-BEC crossover phenomena. 

Fig. Ht^a) shows the crossover phenomena for moderately high densities with (stronger) 
pairing strength J. Here we have set U/t=2. The scaled chemical potential /i' = /i/4t slips 
below the noninteracting band for n = 0.75 and 0.8, in contrast to the existing results where 
BCS-BEC crossover occures for low densities only. 23|] Interestingly, for a higher n (= 0.85), 
the system directly goes from normal to a BEC state without visiting the conventional BCS 
state. 

Next, we explore the BCS-BEC crossover with variation of n for fixed values of U and 
J. For moderately weak parameter values, J/t = 1.5 and U/t = 1, the system prefers to 
stay in the BCS state for the calculated densities (FigJl](b)), consistent with the existing 
results. Interstingly, for the same density values, the system goes to the BEC state as we 
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FIG. 5: (color online) Variation of (scaled) chemical potential, ^' as a function of pairing energy J 
for different system sizes L = 12, 24, and 32 with n = 0.75 and U /t = 2 fixed. Also shown JBEc/t 
with 1/-L, which changes only in the second decimal place as system size is increased. 

increase the pairing strength and Coulomb repulsion, J/t = 4.9 and U /t = 2. Further, the 
system goes to a BCS state quite close to half filling. These results are new and indicate 
the importance of Coulomb repulsion on the crossover mechanism. The reason for choosing 
unusual parameter values {e.g. J/t=1.5 and 4.9) is that, inspite of a thorough and careful 
search for the self-consistent solutions at other parameter values, the success was limited. 
However we could scan a broad density regime which has contributed to a lot more enriched 
study of the problem compared to what already exists. 

To investigate the finite-size scaling effects of the the crossover scenario, we have consid- 
ered lattices of different sizes, namely L = 12, 24, and 32 for fixed value of density n = 0.75 
and Coulomb repulsion U/t = 2. As shown in Fig. [5l the crossover phenomena is qualita- 
tively unaffected as system size is increased. In the inset we have plotted JBEc/t with 1/L, 
where Jbec is the critical value of the pairing strength at which the BCS-BEC crossover 
occurs. It is clear that the crossover phenomena is minimally affected and changes only in 
the second decimal place as system size is increased, which will not affect our qualitative 
study and leave the notion of a BCS-BEC crossover in a d-wave superconductor intact. 

In case of the s-wave superconductors, largely overlapping of Cooper pairs (BCS) 
smoothly evolve into short ranged tightly bound pairs (BEC) with increasing interparticle 
attraction among the fermions. In d-wave systems, the pairs always retain a finite spatial 
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FIG. 6: (color online) Variation of root mean squre radius Rrms as a function of (a) pairing energy 
J and (b) electron density n. Rrms shrinks from a few tens of lattice spacing to order one of lattice 
spacing across BCS-BEC crossover and consistent with the chemical potential picture (Fig. H]). 

extent. A better physical picture of the BCS-BEC crossover can be given by calculating the 
average pair coherence length or the mean pair radius Rrms- It is defined by the following 



relation 



38 



I\ 


mi 


2 ^2^3j„ ^ 




2 


I 




|2 d^r Yl 


1 9k 1 


2 



(19) 



where f{r){g{k)) is the wave function for a Cooper pair in real (momentum) space. 

Fig. shows the evolution of Rrms as a function of the exchange potential and at a 
fixed density in a 24 x 24 lattice. Clearly, Rrms shrinks from a few tens of lattice spacing 
to order one of lattice spacing as the attracting pair strength increases. The transition 
takes place precisely at the same value of J where /i' slips below the noninteracting band 
edge. This supports the proposed crossover scenario in a (i-wave superconductor. For a 
fixed value of exchange potential, the overlap of the electronic states increases with particle 
density, which supresses the emergence of bosonic degrees of freedom near half filling and 
the system instead remains fermionic (BCS-like), even for moderately stronger couplings. 
This is manifested in Fig. El^b), where for n > 0.9 the system shows characteristics of a BCS 
superconductor. In our study, the mean pair radius and the chemical-potential descriptions 
of BCS-BEC crossover are fully consistent with each other (Fig. [Hand Fig. [6]). In contrast 
to the earlier results, where the crossover is seen to occur only in the dilute limit, 23| we 
have seen that it can be realized near half filling. 
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V. CONCLUSIONS 



The BCS-BEC crossover has been investigated for rf-wave superconductors in 2D using 
a t-J-U model near half filling. We have employed the mean-field Bogoliubov-de Gennes 
method and calculated the ground state properties self consistently. To compare and con- 
trast between the zero and large density limits, we included a discussion on the threshold 
exchange potentials that are required to form a two particle bound state (at zero density) 
and superconductivity (at large densities). While the former has no effect on the Coulomb 
repulsion, f/, the latter seems to have been affected rather significantly by it. Further, the 
BCS-BEC crossover scenario is investigated as a function of both the exchange interaction 
and Coulomb repulsion, and the possibility of the crossover is concluded for some specific 
parameter values. To complement the earlier studies on the subject at low densities, the 
crossover picture is demonstrated at high densities by studying the behavior of the chemical 
potential which when falls below the noninteracting band minimum, signals the onset of a 
BEC like phase. A robust support of the crossover is provided by computing the pair radius 
which shrinks from tens of lattice spacing to that of a very few at the onset of crossover. 

It is worthwhile to mention that going beyond the standard mean-field approach by 
including pair fluctuations within the attractive Hubbard model, it was found that around 



half filling the smooth evolution from the BCS to the BEC limits is interrupted. |39l] The 
transition temperature to a superfiuid phase vanishes near half filling over an extended 
range for moderately strong attraction when the system approaches the bosonic regime. 
This vanishing is associated with the divergence of the mass of the pairs or localization of 
the pairs due to pairing fluctuations, which eventually destroys the superfiuid state. Also, 
as the density approaches half filling the chemical potential get pinned to its noninteracting 
value due to particle-hole symmetry. So, although vanishes due to strong interaction, 
the system still stays in the fermionic regime (/i > 0). Thus including fluctuation effects to 
study the crossover phenomena may be a worthwhile exercise for the future. 
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